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ABSTRACT 

We report the discovery of a low-mass companion orbiting the metal-rich, main 
sequence F star TYC 2949-00557-1 during the MARVELS (Multi-object APO Radial 
Velocity Exoplanet Large-area Survey) Pilot Project. The host star has an effective tem- 
perature Teff = 6135 ±40 K, logff = 4.4 ±0.1 and [Fe/H]=0.32 ± 0.01, indicating a mass 
of M = 1.25±0.09 Mq and R = 1.15±0.15 Rq. The companion has an orbital period of 
5.69449 ± 0.00023 days and straddles the hydrogen burning limit with a minimum mass 
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of 64 Mj, and may thus be an example of the rare class of brown dwarfs orbiting at 
distances comparable to those of "Hot Jupiters." We present relative photometry that 
demonstrates the host star is photometrically stable at the few millimagnitude level on 
time scales of hours to years, and rules out transits for a companion of radius ^0.8 Rj 
at the 95% confidence level. Tidal analysis of the system suggests that the star and 
companion are likely in a double synchronous state where both rotational and orbital 
synchronization have been achieved. This is the first low-mass companion detected with 
a multi-object, dispersed, fixed-delay interferometer. 



1. Introduction 



Studies of the frequency, parameter distributions and correlations of extrasolar planets require 
homogeneous samples of hundreds of planets to obtain statistically significant results. Moreover, 
such a sample must have well-understood completeness limits, selection effects and biases, which 
are easiest to obtain from a single, large-scale survey. Given current constraints on the frequency of 
giant planets, detection of such a large sample of planetary systems generally requires a precision 
radial velocity (RV) survey of many thousands of stars. Such a survey also provides a wealth of 
ancillary science. In particular, it is exquisitely sensitive to more massive companions, and because 
it targets a large and broad sample of host stars, it is naturally sensitive to rare binary systems in 
poorly explored regions of parameter space. 

Of particular interest are constraints on the frequency and parameter distributions of low- 
mass companions to solar- type stars with masses near the hydrogen burning limit. One of the 
early results from precise RV searches was the apparent paucity of brown dwarf companions with 
minimum masses (12 Mj < msini < 80 Mj) at separations of a 5-. 5 AU, relative t o more massive 
stellar companions and less massive planetary companions (jMarcy &: Butleiil200d ). Note that we 
denote i as the inclination angle between the companion's orbital angular momentum vector and the 
line-of-sight, and we reserve I as the inclination angle of the stellar rotation axis to the line-of-sight. 
Whil e the frequency of brown dwar f companions at larger separations is still relatively uncertain 
(e.g., iMetchev &: Hillenbrand! l2009l). a meta-analysis of sets of known companions to solar- type 
stars bv lGrether &: Lineweaveii (j2006l ). with corrections for observational bias, confirmed the lack 
of brown dwarfs at close separations. These authors place the 'driest' part of the brown dwarf 
desert at ~ 20 — 50 Mj, with a frequency of companions < 0.5% in this range of masses. 

Although there has b een a steady increase in the number of known brown dwarf candi- 



dates via the RV techniq 
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Marcv et al. 
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Udrv et al. 


20od: 


Endl et al. 
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: Patel et al 
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Wittenmver et al. 


2009: 


Kane et al. 
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Jenkins et al. 


2OO9I: 


Niedzielski et al. 


2OO9I; 


Omiva et al. 



2OO9I ). most of these detections have been at separations a > 0.8 AU. Notable exceptions in- 
clude the transiting brown dwarf CoRoT-Exo-3b with a period of ~ 4 days orbiting an F3V star 
(jPeleuil et al.ll2008l ). and HD41004Bb with a period of ~ 1 day orbiting the M dwarf component 
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of a K-M binary system ([Santos et al.ll2002l ). Brown dwarfs at such short orbital separations are 
of particular interest for several reasons. First, the frequency of such systems as a function of 
their physical and orbital parameters provide diagnostics that may be able to distinguish between 
the various mechanisms that have been in voked for their formation and dynamical evolution (e.g., 
Armitage &: Bonnelll 120021 : iMatzner &: Lev in 2005). In particular, these systems offer observational 
constraints o n the poorly-understood t heory of tidal interactions between host stars and close com- 



panions (e.g.. iMazeh 



2005 



Pont II2OO9I ). Second, these systems are much more likely to transit than 



their longer-period counterparts, as the transit probability is inversely proportional to orbital sep- 
aration. Transiti ng systems yield valuable measurements on th e masses, radii, and mean densities 
of brown dwarfs dstassun et aklbood . boO?! : lOeleuil et ahlboosl ). 



Here we report the discovery of a candidate short-period, brown dwarf companion to the metal- 
rich star TYC 2949-00557-1, a main sequence F star with apparent brightness V ~ 12.1. This 
companion was discovered as part of the Multi-object APO Radial Velocity Exoplanet Large-area 
Survey (MARVELS) Pilot Proje ct (hereafter MP P). The MPP used the W. M. Keck Exoplanet 
Tracker ( Keck ET) instrum ent (jGe et al.l l2006al ) on the Sloan Digital Sky Survey (SDSS) 2.5m 
telescope (jGunn et al.ll2006l ) at the Apache Point Observatory. The Keck ET in strument is a rn ulti 



object (59 targets per exposure), dispe rsed fixed-delay interferometer (DFDI, iGe et al 



2002 



Ge 



2002 



Erskine 



2002 



Erskine et al.ll2003l ). In this instrument, fiber-fed starlight from the telescope 



is first passed through an iodine cell that acts as a stable wavelength reference. This light is then 
fed through a fixed-delay interferometer controlled via a piezoelectric transducer (PZT), and finally 
through a spectrograph that has a spectral resolution of R = 5100. Radial velocity information is 
then imprinted in the phases of the fringes perpendicular to the dispersion axis of the spectrum 
due to a fixed variation in the interferometer delay along this direction. 

The primary goal of the MPP was to demonstrate a fully-integrated DFDI instrument capable 
of observing multiple stars in a single exposure. The second goal was to demonstrate that such an 
instrument is capable of achieving Doppler RV precision sufficient for extrasolar planet detection. 
The final goal was to formulate an operations procedure for conducting a survey using the Keck ET 
i nstrument in a n efficient manner. These three goals were necessary preparations for the MARVELS 
(|Ge et al.ll2009l ) survey: a multi-object, DFDI, extrasolar planet survey that is part of Sloan Digital 
Sky Survey III (SDSS-III)EI 

The MARVELS Pilot Project was conducted in 2007, and consisted of 5-38 observations of 
708 stars taken over 1-5 month baselines. Although the primary purpose of the MPP was to 
lay the groundwork for the full MARVELS survey, the cadence, number, and precision of the 
MPP was nevertheless sufficient to detect massive companions to a nur nber of the target stars. 
The RV data for each star was initially fit using the RVSIM program (IKane et al.l 120071 ) . The 
companion to TYC 2949-00557-1 emerged as an excellent candidate, and here we describe the 
discovery radial velocity data, as well as additional spectroscopic and photometric data acquired 
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to confirm the companion and furtlier characterize the host star. TYC 2949-00557-lb is the first 
low-mass companion detected with a multi-object, dispersed, fixed-delay interferometer; previous 
observations with a single-object DFDI instrument at Kitt Peak National Observatory resulted 
in the first extrasolar planet discovered via this te chnique (|Ge et al.l l2006bl ). as well as the first 
confirmed planet via DFDI (Ivan Evken et al.l 12004 ) and the ability to measure precise, absolute 



radial velocities with DFDI (Mahadevan et al 



20081 ). 



2. Doppler Observations 
2.1. MPP Observations 



The MPP targeted 708 F, G, and K dwarfs with 7.6 < F < 12 in 12 different fields, each 
containing 59 stars. Data for each field was processed simultaneously using a pipeline developed 



for m ulti-object DFDI instruments (see Ivan Evken et al.ll2004l : iGe et al.ll2006bl : iMahadevan et al 
20081 . for details of basic DFDI processing steps.) For each target, we determine a "quality factor'' 
(QF), which we define as 

RMS (X- (X)) 



QF 



(1) 



MEDIAN (dx) ' 

where X represents the RV measurements and RMS is the root-mean-square residual. Targets that 
have QF > 25 either have intrinsically large RV variability, or they are stars with line-of-sight 
rotational velocities >10 km s^^, from which we are unable to extract precise RV measurements. 
TYC 2949-00557-1 was identified from a field with 58 other targets as having a QF of 37.20, in this 
case indicative of its binary nature. A total of 14 usable Doppler RV measurements were obtained 
spanning 116 days from Jan- Apr 2007. Table [Tj contains the barycentric Julian date using the 
Barycentric Dynamical Time standard (BJDtdb); measured Doppler RV and associated errors for 
all 14 observations. The Doppler velocities presented in Table [Tj are absolute velocities calibrated 
to the solar spectrum, and are the sum of the barycentric velocity of the system and the additional 
Doppler variability caused by the companion. In order to provide an indication of the true level of 
systematic errors in the MPP data for objects in this field, we note that the 47 (likely constant) 
targets with QF < 25 in this field have a median quality factor of 3.25. 

We searched for periodicity i n the Doppler data for TYC 2949-00557-1 using a Lomb-Scargle 



(Lomb 


1976: 


Scargle 


1982; 



Gumming et al.lll999l ). The resulting 
power spectrum is shown in Fig. (H revealing a clear peak at P = 5.68 days with a power of ~ 217. 
To assess the significance of this peak, we ran a Monte Garlo simulation with 10^ trials. For each 
trial, we scrambled the times of the data points, computed the periodogram, and recorded the most 
significant peak. We found no trials with power greater than that of actual data, indicating a false 
alarm probability of < 10~^. 



The advantage of a multi-object instrument is that the multiple targets that are observed 
simultaneously can be used to check for common systematics trends in data. Thus, as an additional 
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check on the rehabiUty of the observations, we constructed periodograms for the other 58 objects 
on the plate and calculated the power for each object at the period of the suspected companion. 
Any common systematics present in the data due to the sampling rate or instrumental effects will 
result in significant power at a common period for other targets. None of the other 58 targets have 
significant power at the period of TYC 2949-00557-1; the next strongest candidate has an FAP at 
that period of 77% and the other targets have FAPs > 99%. 

The best-fit amplitude derived from the periodogram is ~ 5500 m s~^, indicating a minimum 
mass in the brown dwarf regime for an FGK-type primary. Given that brown dwarf companions in 
this period range are rare, we decided to obtain additional precise radial velocity measurements, 
high-resolution spectra, absolute photometry, and precise relative photometry time series, in order 
to better characterize the primary and ascertain the nature of the companion. 



2.2. KPNO Doppler Observations 

Observations for the purpose of confirmin g the Doppler v ariability and orbit were conducted 



with the Exoplanet Tracker (ET) instrument (iGe et al 



2006bl ) at the Kitt Peak National Obser- 



vatory with the 2.1m telescope. Two observations separated by several hours were taken each 
night over seven consecutive nights starting on Oct. 10, 2008. A total of eleven usable epochs were 
obtained. Integrations consisted of 60-min exposures bracketed by exposures of a Tungsten lamp 
passing through an iodine gas cell that acts as a calibration for instrument drift. Each arm of 
the interferometer produces a DFDI spectrum from which radial velocities are measured. The two 
beams are processed separately, and their measured radial velocities are combined via a weighted 
average based on the RV uncertainties. Table [2] contains the dates and velocities for the KPNO ET 
measurements. Unlike the results from the MPP, the velocities presented in Table [2] are relative 
RVs, i.e., measured relative to one of the epochs. Note that none of the velocities are exactly zero 
because an instrumental drift has been subtracted off based on the calibration lamps, and they 
are zeroed to a different epoch than the star. Because these data are relative RVs, an offset exists 
between the values in Table [1] and Table [2] that must be included as an additional parameter when 
performing a combined analysis of the two data sets. 



2.3. HET Doppler Observations 



Observati ons of the candidate were also conducted using the R=60,00 mode of the HRS 
spectrograph (jTull Ill998l ) on t he Hobby-Eber l v Tel escope (HET) telescope (jRamsey et al.lll998l ) 
in the queue scheduled mode (jShetrone et al.l 120071 ). The spectra consisted of 46 echelle orders 
recorded on the "blue" CCD detector (407.6-592 nm) and 24 orders on the "red" one (602-783.8 
nm). The spectral data used for RV measurements were extracted from the 17 orders that cover the 
505-592 nm range of the iodine spectrum. A total of ten Doppler RV measurements were obtained 
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spanning 83 days from Dec 2008 through Feb 2009. The starhght was passed through an iodine ceh 
to provide a stable reference to caUbrate instrument drift. Two exposures without the iodine ceh 
were taken to act as stellar templates. Due to the faintness of the target, the RVs were computed 
relative to each template and a mean value was determined. The results for both templates agree to 
within three sigma. Table [3] contains the dates and velocities for the HET measurements. Similar 
to the results from the KPNO ET in Table O these velocities are relative RVs, and therefore an 
offset exists between these values and the ones in both Table [T] and Table [2j 



2.4. Combined RV Analysis 

In order to check for consistency, we first fit the MPP, KPNO and HET datasets individually 
to a seven parameter RV fit, where the seven parameters are the velocity semi- amplitude K, ec- 
centricity e, argument of periastron w, period P, time of inferior conjunction of the companion Tc, 
velocity zero point 7, and linear slope 7 (in order to allow for additional companions or systematic 
drifts). The best-fit solution was found using a hybrid downhill-simplex fit to the nonlinear param- 
eters and an exact (linear) fit to the linear parameters. There are 14 MPP points, 11 KPNO points 
and 10 HET points, so there are seven, four and three degrees of freedom (dof), respectively. 

For the MPP fit, we find a x^/dof of 44. Given that the data points basically follow the 
model, and that the more precise HET data (whose error bars are overestimated, see below) fit the 
model well, the large x^/dof indicates that there are systematic uncertainties in addition to the 
photon noise, and thus the errors are severely underestimated. Given the large median QF=3.25 
found for the majority of the (likely constant) stars in this field, this result is not surprising. 
Indee d, there is a known syst ematic error in DFDI when utilizing an iodine cell in the stellar beam 



path ( van Eyken et al.ll2010l ). For the KPNO fit we find a x^/dof of 6.96, once again indicating 
the uncertainties are underestimated. For the HET fit, we find a x^/dof of 0.09, therefore the 
uncertainties are likely overestimated for the HET dataset. It is worth noting that a statistically 
significant slope is found when fitting the HET data. Fitting the HET data without a slope produces 
a significantly worse x^/dof . 

We performed additional RV fitting to test the significance of the HET slope. Fitting the 
HET data with no slope and eccentricity forced to zero still results in a x^/dof <C 1. Since this is 
the simplest model of an orbiting companion, it confirms the HET uncertainties are overestimated. 
Several different models, in which slope is a free parameter, eccentricity is a free parameter, or both 
are free parameters, all result in lower x^/dof. There is no evidence of nonzero eccentricity in any 
of the best-fit solutions, but a significant slope is found in all cases when left as a free parameter. 
To be self-consistent, we allow for slopes in the MPP and KPNO data fitting as well, and note 
that their best-fit slopes are consistent with the HET value, but are much more poorly constrained 
due to the much larger RV uncertainties in those data sets. As a final check, we fit all data sets 
with eccentricity forced to zero and no slope, and find that the other orbital parameters are not 
qualitatively different from the result where slope is left as a free parameter. We therefore choose 
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the case of zero eccentricity and non-zero slope as our preferred solution. 

Uncertainties in the fitting parameters will be inaccurate if it is determined using misestimated 
RV errors. It is therefore important to attempt to correct the errors such that x^/dof ^ 1. 
However, given that we do not know why the errors are misestimated, particularly in the case 
of underestimated errors, the appropriate method to correct the uncertainties is not clear. Our 
approach was to try several different ways of correcting the errors to force x^/dof = 1. Specifically 
we investigated four different cases for treatment of the MPP (KPNO data come from a similar 
pipeline) and HET RV uncertainties: a scaling of the errors by a constant factor, an addition in 
quadrature of a constant error, a removal of suspected outliers based on the magnitude of the RV 
uncertainty followed by scaling, and a treatment of all data points with a constant error value. 
Ten Markov Chain Monte Carlo (MCMC) simulations were run for each case. The starting values 
for the parameters in these chains were chosen to span a range that is large with respect t o the 
expec ted la uncertainty, and the chains were stopped after reaching convergence as defined in iFord 
(|2006l ). and then the chains were merged. 

Analyzing the MPP and HET data separately yields discrepant periods at the ~ 2a level for 
all cases of error treatment. We also found that the choice of error treatment can affect the derived 
value of ecosw from the MPP data as well at the ~ la level. However, the other parameters from 
the MPP data, as well as all parameter values from the HET data, were not significantly affected 
by different treatments of the uncertainties. From this test, we conclude that there is no strong 
justification for removing any data points from the fit, so we conducted the final joint analysis of 
all three datasets where the each set of errors are scaled by a constant factor. 

The MPP errors are scaled by a factor of 6.64, the KPNO errors are scaled by 2.64, and the HET 
errors are scaled by a factor of 0.3 such that the reduced is ~ 1 when each is fit independently. 
The fit including eccentricity as a free-parameter is consistent with zero eccentricity, therefore we 
run a second fit that forces e = 0. In that case the error scalings are factors of {6.82, 2.71, 0.24} 
for MPP, KPNO and HET, respectively. We further scale the errors of all three data sets by a 
factor of 1.30 for the case where eccentricity is left as a free parameter, and 1.25 for the case 
where eccentricity is fixed to zero. This scaling is done so that the x^/dof = 1 in the combined 
fit, and is necessary due to the systematics present in the data sets. Given the close separation 
of the companion, it is expected that the orbit has been tidally circularized, consistent with our 
findings. We therefore treat the case with eccentricity fixed at zero as our final values, but quote 
the parameters from both cases in Table HI which contains the values of the orbital parameters for 
the case of non-zero eccentricity (eccentric) and eccentricity fixed at (circular). 

Fig. [2] shows the final results of the joint RV fitting and fixing the eccentricity at 0. MPP data 
are the blue squares, KPNO data are the green triangles, HET data are the red circles, and the 
systemic velocity 70 has been removed. We find 70 = 18.68 it 0.24 km s~^ for the star, with an offset 
between the MPP and KPNO data of 14.90 ± 0.25 km s"^ and an offset between the MPP and 
HET data of 12.61 ± 0.24 km s"^ The final orbital period is determined to be 5.69449 ± 0.00023 
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days and an RV semi-amplitude of 6.113 it 0.009 km s~ . We searched for an additional signal in 
the residuals from the joint fit that might be caused by an additional companion in the system, 
but found no other frequencies with significant power. 

3. Relative Time Series Photometry of the Host Star 

Photometric observations are an important step in analyzing low-mass companions discovered 
via the Doppler technique. High-precision photometry can be used to search for transits of the 
companion. Additionally, time-series photometry can be used to rule out stellar mechanisms of 
Doppler variability, such as chromospheric activity due to starspots or stellar pulsations. In the 
case of stars with detectable starspots, time-series photometry can be used to determine a stellar 
rotation rate. In this section we present and analyze time series relative photometry of TYC-2949- 
00557-1 from two sources: relatively precise (few mmags) photometry covering a relatively short 
timespan (2-8 hours over five nights) from the Hereford Arizona Observatory, and less precise (few 
percent), but more comprehensive photometry consisting of 7194 epochs taken over roughly three 
years as part of the Kilodegree Extremely Little Telescope (KELT) North transit survey. Neither 
datasets show any evidence for variability of the host star. 

3.1. Relative Photometry from Hereford Arizona Observatory 

Initial photometric observations of the primary were performed on four nights in 2009 (2/19, 
2/21, 2/27, and 3/16) at the Hereford Arizona Observatory (observatory code G95 in the lAU 
Minor Planet Center) , a private facility in southern Arizona. Additional observations were made in 
2010 on 4/15 to search for transits based on an updated transit ephemeris from the combined RV 
analysis ( ^2.4p . All data were taken with an 11-inch Celestron Schmidt- Cassegrain (model CPC 
1100) telescope that is fork-mounted on an equatorial wedge, an SBIG ST-8XE CCD with a KAF 
1602E detector, and an SBIG AO-7 tip-tilt image stabilizer used to maintain the field at a fixed 
position on the CCD. The observations in 2009 were done without a filter ("C" band), resulting in 
an effective central wavelength of ~570 nm between Johnson V and R bands. The observations in 
2010 were taken with a Sloan r' filter. Data toward the end of the night on 4/15/2010 were taken 
at very high airmass (out to sec 2: = 5.7), resulting is somewhat degraded photometric precision. 

Fig. [3] shows the relative photometry over the five nights, which demonstrate that the primary 
star is intrinsically stable on the time scale of several hours, at the level of 2-4 millimagnitudes. 
Based on the final ephemerides determined in §2.41 only the last night (4/15/10, top row) covers 
possible times of predicted transits. The vertical bars are the predicted times of ingress, mid-transit 
and egress based on the RV fit in §2.41 for the two methods of RV fitting ("C" is for e = 0, "E" is 
for non-zero eccentricity), and an assumed transit duration of 3.3 hours, corresponding to a nearly 
central transit. The widths correspond to the uncertainties in the mid-transit times. There is no 
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evidence of a transit at the most likely depth (~ 0.8%) and duration during these observations. 
In ^3.31 we consider the uncertainties in the ephemeris and properties of the primary, as well as a 
range of impact parameters, to quantify the confidence with which we can exclude transits in this 
system. 



3.2. Relative Photometry with KELT 



The Kilodegree Extremely Little Telescope (KELT) North survey is a wide-field photometric 
survey of ~ 40% of the northern sky designed to monitor fairly bright (8 < V < 12) stars in order 



to search for planetary transits ( Pepper et al. 2003 . 2007 : Siverd et al 20091 ). The KELT survey 



instrument consists of an Apogee AP16E (4K x 4K 9//m pixels) CCD camera attached to a Mamiya 
645 medium-format 42mm aperture camera lens. The resultant field of view is 26° x 26° at roughly 
23" / pixel. The standard configuration uses a Kodak Wratten ^8 red-pass filter and the resultant 
bandpass resembles a widened Johnson-Cousins i?-band. KELT-N is permanently mounted on a 
fixed pier at Winer Observatory in Sonoita, AZ. 

The KELT-N survey targets 13 star fields centered at 31.7° declination (the survey site latitude) 
spaced fairly evenly through all 24 hours of R.A. with slight overlap. Exposure times are 150 
seconds, which yields relative photometric precisions of better than a few percent for V < 12. 
Typical cadences are roughly 20 to 30 minutes when the target field is visible, and to date there 
exist ~ 5000 — 7000 epochs per target. The areal sky coverage, target magnitude range, and high 
photometric precision of the KELT survey results in excellent synergy with the MPP (as well as 
the full MARVELS survey). TYC 2949-00557-1 is in one of KELT's target fields and is a good 
example of this synergy. We use the KELT photometry to characterize the photometric variability 
of the host star, and to search for signatures of transits of the companion. We first briefly describe 
the data reduction, and then describe the light curve analysis. 



Images of the field ar e fiat-fielded, and then relative photometry is extracted using the ISIS 



imae 



,ge subtraction package ( Alard Sz Luptonlll998 ). in combination with DAOPHOT (Stetson 1987 ) 



to perform point-spread function fitting photometry. We further eliminate problematic images due 
to poor observing conditions by examining outliers from the ensemble of individual light curves 
on the CCD. Any epochs that produce photometric outliers in a significant frac tion (>, 5%) of 



light curves are removed from all light curves. We use the VARTOOLS program (jHartman et al 



20081 ) to remove co mmon trends due to systematic errors from the light curves using the Trend 
Filtering Algorithm (jKovacs et al.ll2005l ). first removing the 20 points with the highest flux and 
the 20 points with the lowest flux from all light curves. We choose the 400 stars with the lowest 
RMS values as comparison stars for trend removal, ensuring that these stars are evenly distributed 
across the region of the CCD near the target, and excluding variable stars and saturated stars. 
Finally, the errors of all the light curves are scaled by a constant value, chosen to force the modal 
value of x^/dof for a weighted, constant flux fit to the light curves to be unity for the ensemble. 
For TYC 2949-00557-1, this procedure resulted in a x^/dof that differed somewhat (< 50%) from 
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unity. Since we see no evidence for photometric variability for this star (see below) we further scale 
the errors to force x^/dof = 1 to be conservative. TYC 2949-00557-1 happened to fall in an overlap 
region of two target fields, and as a result we had two sets of light curves for the target. These 
data were reduced independently, and then combined after subtracting the difference between their 
weighted mean magnitudes. 

The KELT lightcurve for TYC 2949-00557-1 is shown in Fig. SI It contains 7194 data points 
spanning 3.2 years, and has a weighted RMS of 3.5%. We search for variability using a weighte d 



Lomb-Scargle periodogram with floating mean (jLombl 119761 : IScargld Il982l : ICumming et al.l Il999l ) , 
and find no significant peaks (see Fig. and in particular no evidence for periodic variability 
near the period of the companion {P ~ 5.69 days), or the first harmonic {P/2). Fig. [5] shows the 
KELT lightcurve, phased according to the best-fit RV ephemeris (Table [4]), as well as binned 0.025 
in phase (3.46 hrs). The RMS of the binned data is ~ 2.4 mmag, and the x^/dof = 0.91 for a 
constant fit, indicating a low level of correlated noise, and no evidence for variability at the few 
mmag level. We limit the amplitude of any variability at P/2 to be < 2 mmag; unfortunately this 
i s well above the level of ellipsoidal variability expected for this companion of ~ 0.03 sini mmag 
dPfahl et allbood ). 



3.3. Limits on Transits 

Given the relatively high a priori transit probability of TYC 2949-00557-lb of ^ R/a ~ 8%, 
where R is the stellar radius, we searched for transits in the KELT dataset combined with the 
Hereford data from 4/15/10. The expected transit duration is ~ RP/{7ra) ~ 3.4 hrs for a central 
transit, and the expected fractional depth is (5 ~ (r/R)'^ ~ 0.8% (r/Rj)'^, where r is the radius of 
the companion. The expected radius of the companion depend s on its true mass, as well as the 



age of the system, but is likely to be ~ 1 Rj (iBaraffe et al.ll2003l ). Unfortunately, while the KELT 



data has excellent phase coverage, it is not of sufficient quality to detect or rule out the expected 
signal. If transits were present, we should expect to detect them with a signal-to-noise ratio of 

1 /2 

S/N ~ ivV2 (Ji) i (2) 
yvra/ a 

where A'^ = 7194 is the number of data points, and a ~ 3.5% is the typical uncertainty. Thus, 
S/N ~ 3(r/i?j)^, which is marginal unless the companion has a radius significantly larger than 
Jupiter. On the other hand, the majority of the Hereford data is generally of sufficient quality 
to detect transits at the expected depth, and one night covers the predicted transit time for the 
companion. Unfortunately, there is no indication of a transit at the expected time. 

We nevertheless proceed with a quantitative search for a transit signal. We combine the 
KELT data with the Hereford data from 4/15/10 after first subtracting the difference between 
their weighted mean magnitudes. The slight difference in passband between the HAO and KELT 
data do not affect the ability to detect a transit in the combined data set. We use the distribution 
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of companion periods P and expected transit times Tc from the MCMC analysis of the combined 
radial velocity data described in ^2.4[ For each combination of Tc and P (i.e., for each link in 
the Markov Chain), we draw a random value for the T^s, log 5, and [Fe/H] of the primary from a 
Gaussian distribution, with the central values and dispersioi is given in Tabled as determined from 
the spectroscopic analysis described in ^5.2[ We then use the lTorres et al.l ([20101) empirical relations 
to estimate the mass M and radius R of the primary for those values. We add an additional offset 
to M and R drawn from Gaussia i is wit h dispersions equal to the dispersions of the fits to the 
empirical relations in [Torres et al.l (|2010l ): specifically 6.4% in M and 3.2% in R. Finally, we draw 
a random value of the impact parameter of the transit in units of the radius of the star in the 
range [0,1]. Assumi ng a radius for the com panion, we then compute the expecte d transi t curv e 
using the routines of lMandel &: Agoll (j2002l ). using limb-darkening coefficients from IClaretl (j2000l ). 
assuming that both the KELT and Hereford bandpasses roughly correspond to R. We then fit this 
curve to the combined dataset, and compute the improvement in relative to a constant flux fit 
to the data. We repeat this for each link in the Markov chain, as well as for a variety of different 
companion radii. 

We search for significant improvements in which would be indicative of a detection. Our 
best-fit has Ax^ = —11.7 relative to a constant flux fit. In order to asses the significance of this 
improvement in x^i we repeat the search for "anti-t ransit", i.e., signa ls with the same shape as a 



transit but corresponding to an increase in flux (see iBurke et al.l |2006| ) , and find improvements in 
at similar levels. We therefore conclude that there is no evidence for a transit in the combined 
KELT and Hereford data. 

Given that we have not detected any evidence of transits, we now ask what the probability 
is that we would have detected a transit of a given radius, assuming that the companion transits 
(i.e., b < 1). To do this, we simply determine what fraction of the steps in the Markov Chain 
described above result in an increase in above a certain level, as a function of the radius of the 
companion. This result is shown in Fig. [6l for Ax'^ = {9, 16, 25}. The Ax^ values were chosen as 
representative values: Ax^ = 16 is the likely detection limit, 25 is chosen as a conservative limit 
and 9 is chosen to straddle the true detection limit. The black, long-dashed line is a case where 
flat-bottomed, boxcar-shaped transits (no ingress/egress) were used and represents Ax^ = 16. It 
shows that detailed modeling of limb darkening and the ingress/egress has little effect on the final 
results of this test. Given the properties of the noise as estimated from the improvements in the 
fits from "antitransits" , signals with A^^ ^ 16 are likely to have been reliably detec ted. Thus, we 



can ex clude transits of a companion with r > 0.75 Rj. at the 95% confidence level. iBaraffe et al 
(|2003l ) predict radii of > 0.77 Rj for brown dwarfs with m ~ 60 Mj for ages of < 5 Gyr. We 



conclude it is unlikely this companion transits, unless the system is substantially older than 5 Gyr, 
which is unlikely given the effective temperature and surface gravity of the host star (^. 
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4. Absolute Photometry 



The Tycho-2 catalog's (|H0g et al.l 120001 ) F-band magnitude for this object is 11.840, however, 
Tycho magnitudes are known to significantly degrade beyond > 11. Measurements were taken 
at Hereford Arizona Observatory using both B,V, Rc, Ic and Sloan g',r',i' filter sets. A total of 
64 Landolt standard stars and 15 SDSS standard stars were used as calibrators. The conversion 



equations of lSmith et al.l ([2002) were used to convert the Sloan filter measurements into the BVRI 
system. The agreement between the observed BVRI magnitudes and those converted from Sloan 
magnitudes agree within the measurement error. We adopt the unweighted average of each BVRI 
measurement and the larger of the two statistical errors for the final magnitude results. Table [6] 
summarizes the measured magnitu des and final results for the multi-band photometry. Near-IR 
fluxes are taken from the 2MASS (jSkrutskie et al.l 120061 ) Point Source Catalog and are presented 
in Table O 



5. Characterization of the Host Star 
5.1. SED Analysis 

We use the BVRI fluxes along with the 2MASS near-IR data to flt mode l spectral energy 



diagrams (SEDs) to derive approximate stellar parameters. NextGen models from lHauschildt et al 



(|l999l ) in grids of 100 K for T^s, 0.5 dex for log (7 and 0.5 dex for the metallicity, represented by 
[Z/H], are fit along with the line-of-sight extinction. There is a degeneracy between line-of-sight 
extinction and derived T^s when fitting SEDs with unknown extinction. Fig. [7| (top) shows the 
map in Teg — Ay space. The interior contours represent 1-cr uncertainties assuming that Tgff and 
Ay are the only two free parameters. The exterior contours are the 1-a uncertainties with all four 
parameters {T^g, Ay, logg and [Z/H]) as free-parameters. We find the global minimum in x^-space 
is for a T^e = 6000 K, logg of 5.0, [Z/H] of 0.0 and extinction Ay = 0.5. Fig. [7] (bottom) shows 
the best-fit NextGen model along with the observed fiuxes. Other solutions exist within the 1-a 



contours at cooler Teg and smaller Ay, however, the dust maps of ISchlegel et al.l () 19981 ) give an 
Ay of 0.45 along this line-of-sight, consistent with the best-fit Ay of 0.5 assuming that the star is 
behind the majority of the dust along this line of sight. The hotter temperature solution is also 



consistent with the Tefj derived from Echelle spectra presented in §5.2 



5.2. Spectral Synthesis 



In order to derive physical properties of the host star and estimate the minimum mass of 
the companion, stellar templates from the HET observations that do not contain iodine lines 
were used to derive pa rameters of the host star. We use the latest MARCS model atmospheres 
(jGustafsson et al.ll2008l ) for the analysis. Generation of synthetic spectra and the line analysis were 
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performed using the turbospect rum code (lAlvarez fc Plea Il998l ). which employs line broadening 
according to the prescription of iBarklem &: Q'Maral (|l998l ). The line lists used are drawn from a 
variety of sources. Updated atomic lin es are taken mainly from the Vienna Atomic Line Database 
(VALD) database dKupka et al.lll999l^. The molecular species CH, CN, OH, CaH, and TiO are 
provided by B. Plez (see lPlez Sz Cohenll2005l ). while the NH, M gH, and C2 mol ecules are from the 



Kurucz linelists. The solar abundances used here are the same as 



Asplundl torn . The Fel and Fell 



lines used for the line ana lysis were comp iled by Santos et alJ ( 20041 ) , who use solar Fe abundance 
to derive gf values. Both lAsplundl (j2005l ) and ISantos et al.l (|200J) use HARPS solar spectra and 
an Fe abundance of 7.45. 

We derive a Toff = 6135 it 40 K and [Fe/H] = 0.32 it 0.01, based on Fel excitation equilibrium 
and a logg' = 4.4 it 0.1 based on the ionization equilibrium of Fel and Fell lines and by fitting the 
wings of the Mgb triplet at 5167, 5172 and 5183 A. The error estimates are based on the equivalent 
width of Fe lines and the errors of Fe abundances from the individual lines. A microturbulence 
value = 1-65 km s^^ is derived by forcing weak and strong Fel lines to give the same abundances. 
Fitting the Fe lines in the Mgb region yields a rotational velocity of vsinl = 7 km s~^. We only 
used the 110 Fel lines weaker than 100 mA for the analysis. Fig. [8]shows the continuum normalized 
spectra in black and the best-fit model in red for the Mgb region. 

In addition to the HET t emplates, we obtai ned high-resolution (R ~ 31,000) Echelle spectra 
using the ARCES instrument (jWang et al. 1120031 ) on the APO 3.5m telescope. Seven exposures for 
a total of 63 minutes of integration were obtained. Data was reduced using a modified IRAF script 
originally written by J. Barentine and J. Krzesinski for ARCES data. Spectra are corrected for bias 
and dark subtraction, cosmic rays and bad pixels. Flatfielding is performed using a quartz lamp 
and two different sets of integration times: a "blue" set of 4-min integrations using a blue filter and 
a "red" set of 7-sec integrations with no filter in the beam. These sets are then combined to form 
a master flatfield image. The two different quartz sets are used to maximize the signal-to-noise in 
both the blue and red end of the spectrum. Spectra are wavelength calibrated using a sequence of 
10-sec ThAr integrations taken a few times during each night. The star HD 42088 (spectral type 
06V) was observed to remove telluric lines. 

The spectrum was ana lyzed using the IDL-based program Spectroscopy Made Easy, or SME 
(jValenti &: Piskunov Ill996l ). This code uses synthetic spectra and multidimensional least-squares 
minimization to determine the best set of stellar parameters for an observed spectrum. These pa- 
rameters include effective temperature, surface gravity, metallicity, microturbulence, macroturbu- 
lence, projected rotational velocity, and the radial velocit y. We follow the guidelin e s from previous 



spect roscopic studies of host stars that used SME (e. s.. IValenti &: Fischer I l2005l : IStempels et al.. 
20071 ). We used three-dimensional interpolation on the lKuruczl(ll993l l grid of local thermodynamic 



equilibrium (LTE) model atmospheres and the VALD database to obtain line data for transitions 
with predicted absorption cores deeper than 0.5% of the continuum. For the VALD queries, we 
used solar abundances, Tcs = 5770 K, logg = 4.44 and = 0.866 km s~^. 
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For the SME analysis, as suggested by IStempels et alj (j2007l ). we fixed the parameter to 



0.85 km s ^, in order to decouple the correlation between n iicroturbulence (.t and me tallicity. For 



the macroturbulence (t, we follow the empirical relation of IValenti Sz Fischer I (j2005l ) which gives 
Q = 4.5 km for a star with Teg ~ 6200 K. We were unable to obtain consistent results from 
the gravity-sensitive Mgb triplet region, therefore we fix logg' at three values of 4.3, 4.4, and 4.5, 
corresponding to the l-a range determined from the HET spectra. We set as free parameters Tefj, 
[M/H], f sin/, and the radial velocity Vrad and utilize the metal-rich region of 6000 — 6200 A. The 
uncertainties of the parameters are derived from the range of best-fit results using the three fixed 
log 5 values. We derive T^s = 6246^45 K, [M/H] = 0.3615+°;°°^ and wsini = 7 ± 1 km s'^, in 
reasonable agreement with the values derived from the HET spectra. Fig. [9] shows the best-fit 
model in black and the input spectrum in white for a portion of the 6000 — 6200 A range used in 
the fitting. 

The SME analysis of the ARCES spectrum is used as an independent check on the derived 
temperature and metallicity from the HET spectra. Because it is based on a smaller wavelength 
region, and cannot be used to independently derive the surface gravity, we chose to adopt the 
results from the HET analysis for the final stellar properties. 



6. Determination of Host Star Mass and Radius 



An alternative to interpol ating isochrone mod els to determine stellar properties is to apply the 
analytical equations derived bv lTorres et al.l (|2010l ) using measurements of eclipsing binary systems. 
We use this empirical relation to derive the mass and radius of the primary using the values of T^s, 
log 5 and [Fe/H] obtained from the HET template spectra. Applying the equations for logM and 
logi? yields a mass of M = 1.25 zfc 0.09M(t^ an d a radius of R = 1.15 it O.IBRq. Correlations of 



the best-fit coefficients from lTorres et al.l (|2010l ) are included in the errors, but correl ations of Tpff, 
log Q a nd [Fe/H] are not considered. The reported scatter in the relation as found in lTorres et al.l 
(|2010l ) of (Tiogm = 0.027 and aiogr = 0.014 are also included in the mass and radius uncertainties, 
respectively, by adding them in quadrature. 

Fig. [TOl compares the spectroscopically measured T^s and log 5 of TYC 2949-00557-1 (red 
erro r bars) to a theoretica l stellar evolutionary track from the Yonsei-Yale ("Y^") model grid 
(see iDemarque et al.l l2004l . and references therein) . The solid curve represents the evolution of 
a single star of mass 1.25 Mq and metallicity of [Fe/H]=+0.32, starting from the zero-age main 
sequence (lower left corner), across the Hertzsprung gap, and to the base of the red- giant branch. 
Symbols indicate various ages along the track labeled in Gyr. The dashed curves represent the same 
evo lutionary track but for masses ±0.09 Mq, representative of the la uncertainty in the mass from 
the [Torres et al.l (|20ld ) relation. The filled gray region between the two mass tracks represents the 
range of expected locations of a star like TYC 2949-00557-1 given the 1-a mass uncertainty and 
measured metallicity. We emphasize that we have not dir ectly measur e d the mass of TYC 2949 
but rather derived the mass using the empirical relation of [Torres et al.l (I2OIOI ). Our purpose here 
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is not to test the accuracy of the theoretical stellar evolutionary tracks, but rather to constrain 
the evolutionary status of the TYC 2949 system. The spectroscopically measured T^g, log (7, and 
[Fe/H] place TYC 2949 near the zero-age main sequence, with an age of at most ~2 Gyr. 

The distance to the host star can be computed once the bolometric luminosity is known. We 
use the Stefan-Boltzmann l aw to derive the lum inosity of the star using the T^g found in §5.21 and 
the radius calculated by the [Torres et al.l (|20ld ) relation. The absolute V magnitude is then given 

by 

Mv = -2.5 log ^ + MboI^, - BCv (3) 



-2.h\og — + MboIq 

where MboIq is the bolometric absolute magnitude of the Sun, here assumed to be 4.74, and BCy 
is the bolome tric correction to the V ba nd. We adopt a value of BCy = —0.175, interpolated from 
Table 15.7 in brihing k Landoltl ([2OOJ) based on the T^s of TYC 2949-00557-1. The distance can 
then be calculated via the distance modulus and assuming a valu e for the line-of-sight extinction. 
If we assume Ay = 0.45 it 0.1, which is consistent with both the ISchlegel et al.l (jl998l ) dust maps 
and the best-fit value obtained from the SED fitting in §5.1^ we derive a distance to the system of 
413 pc. This distance is consistent with our implicit assumption that the star is behind the 
majority of the dust along this line of sight, for likely values of the thickness of the dust layer. 
The quoted uncertainty in the distance includes the uncertainties in the stellar radius, effective 
temperature, line-of-sight extinction and apparent V band magnitude. 



7. Companion Mass 

Given the orbital parameters from the joint radial velocity fit ( ^2.4p . and the estimate of the 
mass of the host star derived from the spectroscopic parameters ( §5.2p . we can estimate a minimum 
mass for the companion of 

mmi„ = 64.3±3.0Mj, (4) 

where we have assumed a circular orbit. This minimum mass is based only on the RV and stellar 
parameters, and ignores the fact that edge-on configurations are likely excluded given the lack of 
transit signature. Since the minimum mass is below the hydrogen burning limit (~ 80 Mj), one 
might be tempted to categorize this object as a brown dwarf. However, since the orbital inclination 
i is unknown, we do not know the true mass, and given that the minimum mass is not far below the 
hydrogen burning limit, it does not suffice to make the usual assumption that the minimum mass 
of the object can be used to characterize its nature. Rather, we must be somewhat more careful to 
estimate the probability distribution for its true mass. 

We proceed to estimate the probability distribution of the companion mass using a similar 
method as was used to search for and exclude transits in the KELT photometric data as described 
in ^ We use the distribution of companion periods P, and velocity semi-amplitudes K from the 
Markov chain derived from the MCMC analysis of the combined radial velocity data described in 
§2.4( for the fit assuming zero eccentricity. For each set of these parameters (each link in the chain) , 
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we draw a random value for the spectroscopically determined primary parameters Tgg, logg, and 
[Fe/H], according to a Gaussian distribution centered on the best-fit values and with dispersions 
equal to the uncertainties (given in Table [5|). We use these v alues of Tpff, l ogg, and [Fe/H] to 
estimate the primary mass M using the empirical relation of [Torres et al.l (j2010l ) . We add an 
additional offset to M drawn from a Gaussian with dispersion equal to the dispersion in the fit to 
this empirical relation (6.4%). We draw a random value of cosi from a uniform distribution in the 
range (0, 1), and then solve for the mass m of the secondary (note we do not assume that m <^ M). 

We weight the resulting distribution of m by a prior on the luminosity ratio I and a prior on the 
mass ratio q. Specifically, we assume luminosity ratios of / > 0.1 are excluded by the lack of features 
due to the companion in the high-resolution spectra. We assume a flux ratio relationship of the form 
I oc q'^'^, as is roughly appropr iate for main-sequen ce stars. The exponent of 4.5 is derived by fitting 
the values found in Table 1 of lTorres et al.l (|201Cll ) for stars with 0.5 < Mq < 1.5. We note that the 
precise value of the exponent for this relationship does not significantly affect our conclusions for 
reasonable values in the range 2.5-6.5. We therefore assume the companion is on the main sequence 
and is not a remnant. Since the mass ratio distribution for companions in the relevant range of 
masses is uncertain, we adopt three different priors tha t are likely to bracket the tru e distribution 



(see, e.g.. iDuquennov Mavor 



1991 



Mazeh et al. 



1992 



Grether &: Lineweaverll2006l ). Specifically. 



we assume linear (diV/ d log q oc q), logarithmic (diV/ d log q oc log q) , and constant {dN/ d log q = 
constant) priors. For the logarithmic and constant priors, a brown dwarf companion is slightly 
favored (~ 66% and 61%, respectively), whereas for the linear prior, a stellar companion is slightly 
favored (57%). Fig. [TT] plots the cumulative probability that the mass of the companion is less than 
a given mass for the three different assumed priors. 



8. Effects of Tides 



Given the proximity of the companion to the host star, tidal interactions could be important in 
this system. The tidal effects between solar mass stars and very low mass stars/brown dwarfs are not 
well-studied. Moreover, tidal models themselves are complicated and uncertain, and observational 
constraints are few, especially among brown dwarfs and very late M dwarfs. The situation is further 
complicated by the ambiguity of the secondary's mass, which partly determines the tidal evolution. 
Nonetheless, we consider the tidal evolution of the system in this section. 



Following iGoldreich &: Soted (119661 ) . we assume that the tidal response of a body can be 
adequately modeled by a single parameter, Q' , which is related to the angle between the lines 
con necting the centers of the two bodies, and the center of the deformed body to its tidal bulge. 
See lHeller et al.l (|20ld ) or lFerraz-Mello et al.l (|2008l ) for recent reviews of tidal theory. Furthermore, 
we require that as the orbit of the secondary evolves, the parameter Q' remains constant. If we 
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assume the orbit is circular, it may be shown that the semi-major axis of the secondary decays as 




^11/2 



(5) 



(jGoldreich Soteiill966l ). where Q'^ is the tidal quality factor of the star divided by two-thirds of 
its Love number, M is the mass of the primary, R is the radius of the primary and m is the mass 
of the secondary. Subscript "*" refers to the primary. 

Equation [5] only applies if the primary's rotation period is longer than the orbital period. 
Should the two be equal, the tidal bulges align, the torques reduce to zero, and the orbital evolution 
effectively halts. There still may be evolution due to the obliquity tide raised on the star, but that 
evolution is orders of magnitude slower, and we ignore it here. As the secondary decays, its angular 
momentum, Lgj-fe = mV GMa, decreases and is passed to the star's rotational angular momentum, 
= C^^MI^Q., where Q. is the prir nary's rotation fr equency, and C* is the primary's moment of 



inertia coefficient. We set C* to 0.1 (jMassarottil l2008l ) . Therefore, over a given time, the change in 
the primary's rotational frequency is 



AL 



orb 



C.MRV 



(6) 



where ALo^b is the change in orbital angular momentum over the same time. Eqs. [5] and [6] can 
be solved together to determine when the rotational frequency equals the orbital frequency (i.e. 
synchronization) when orbital decay stops. Here we ignore the case of the primary's rotation period 
being shorter than the orbital period, but in that case, the secondary will spin down the star, while 
its orbit expands, also driving the system toward synchronization. 

We have examined the tidal evolution of t his system in the range 10^ < Qj < lO^''. Most 



19941 : 



Jackson et al 



2008. but see also 



studies find Qj values i n the range 10^ - 10^ (iMathied I 
Barker fc Ogilvie I ()201Cl )). but in reality this parameter is very poorly constrained. Our chosen 
range covers all plausible values. We also vary i, the inclination of the secondary's orbit to the 
line of sight, for 90° - 1° (from edge-on to nearly face-on) and adjust the mass accordingly. We 
determine the primary's rotation period via the measured value of v sin / and the derived radius of 
the primary, R. For this analysis, we set the primary's equator in the same plane as the secondary's 
orbit (i.e., we assume i = I), but this decision does not qualitatively affect our results. 

In Fig. [12] we show the synchronization times from Eqs. [SHSlover the parameter space defined 
above. We consider two models, the best-fit set of parameters (solid contours), and one in which 
f;sin/ = 4 km s~^ (the 3-sigma minimum value), M = 1.22 Mq, and R = 1.02 Rq (dotted 
contours), i.e. a pathological case, permitted by the observational uncertainties, with values of 
vsinl, M, and R that minimize the tidal evolution. First, note the convergence of contour lines 
at i = 43° (solid contours) and i = 25° (dotted contours). These singularities occur because i 
is small enough that the system is currently synchronized, i.e. vsinl equals the circular velocity 
of the companion. For i greater than these values the secondary is spinning up the primary. We 
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find that for a wide range of and i combinations, the secondary quickly spins the primary up 
to synchronization. For the best-fit and Q'^ ^10^) the synchronization time is ^0.1 Gyr. The 
photometry shows no significant variabihty due to starspots that would indicate a particularly 
young F star, so it is likely that the age of the system is larger than the synchronization timescale. 
We interpret this short timescale as evidence that the secondary has already spun the primary up 
to synchronization. For the pathological case, the time to synchronize could be a factor of a few 
larger, but still small compared to the likely age of the system. We do not show the analogous case 
for large v sin /, M and R, as that case is already synchronized for all values of i. 

We conclude that this system is most likely in a double synchronous state in which the orbit is 
circular, and both primary and secondary rotational frequencies are equal to the orbital frequency. 
We emphasize that this conclusion is tentative; should more information regarding or the relative 
orientations of the primary's spin axis and orbital plane be determined, our analysis will need to 
be updated. 

9. Conclusion 

We have discovered a brown dwarf candidate around the main sequence F star TYC 2949- 
00557-1 during the MARVELS Pilot Project, a wide-area, multi-object, radial velocity search for 
planets. We have characterized the properties of the host star and the dynamics of the orbiting 
companion. The companion straddles the hydrogen burning minimum mass with an m„i[^^i = 
64 Mj, and the orbital period of ~ 5.7 days places this candidate in the brown dwarf "orbital 
separation desert". This desert extends to very low-mass stars (VLMS) with masses up to ~ 150 
Mj, reflective of the fact that the formation and dynamical evolution of these objects is independent 
of the precise substellar mass limit. Therefore, TYC 2949-00557- lb is still a "desert dweller" for all 
but the most extreme face-on inclinations. Tidal analysis suggests that the system is in a double 
synchronous state where both companions have achieved rotational synchronization in addition to 
achieving orbital synchronization. 

If the companion's true mass is substellar, this is an example of a rare class of brown dwarfs 
at orbital distances comparable to those of "Hot Jupiters". Our photometric observations show 
no evidence of variability. Although we are able to exclude transits for likely companion radii, 
the orbital inclination is otherwise unconstrained, and thus there remains the possibility that this 
system is a VLMS companion with a nearly face-on inclination. Depending on the assumed priors 
for the distribution of binary mass ratios, we estimate a probability of 43-66% that the companion 
is below the hydrogen burning minimum mass. The large number of stars surveyed by multi- 
object, radial velocity instruments may allow for additional rare systems to be discovered, including 
brown dwarfs in the mass or orbital separation deserts. Such rare systems can provide important 
constraints to models of the formation and internal structure of low-mass stars and brown dwarfs, 
in addition to elucidating the formation and evolution of extrasolar planetary systems. 
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Table 1. MPP RV Observations. 



BJDtdb 


RV 
(m s-l) 


(m s-l) 


2454101.69079 


13339 


94 


2454105.75520 


19819 


98 


2454106.70228 


14981 


96 


2454128.62211 


19018 


143 


2454128.86491 


17170 


91 


2454130.83623 


13299 


94 


2454136.64529 


15043 


85 


2454136.85894 


15854 


139 


2454163.72245 


14882 


91 


2454188.69749 


21003 


136 


2454191.68122 


15421 


97 


2454194.69216 


21197 


115 


2454195.68446 


24453 


102 


2454217.61115 


22758 


79 



"Errors are not scaled to account for 
systematics. 



Table 2. KPNO/ET RV Observations. 



BJDtdb 


RV 






(m s-i) 


(m s-l) 



2454749.97070 


266 


58 


2454751.88336 


1316 


55 


2454751.94579 


1985 


50 


2454752.89657 


7580 


54 


2454752.95821 


8026 


58 


2454753.90161 


9651 


77 


2454753.96355 


9548 


65 


2454754.90388 


4739 


53 


2454754.96580 


4502 


52 


2454755.91284 


-694 


47 


2454755.97519 


-996 


52 



^Errors arc not scaled to account for 

systematics. 
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Table 3. HET/HRS RV Observations. 



BJDtdb 


RV 

(m s-l) 


(m s-i) 


2454807.74076 





69 


2454808.98051 


4674 


67 


2454825.69507 


2469 


107 


2454829.70100 


2482 


66 


2454881.75734 


-23 


286 


2454882.78369 


3178 


51 


2454883.76208 


9451 


63 


2454887.75592 


165 


48 


2454889.74685 


10881 


51 


2454890.75523 


11513 


47 



''Errors are not scaled to account for 
systematics. 



Table 4. Best-fit dynamical properties of TYC 2949-00557-1. 





\aluc 


1 "iK*(*ri aiiil\" 


\alu(- 


1 "iK-(>rl aiiil\" 




ECCENTRIC CASE 


CIRCULAR CASE 


Period (days) 


5.69459 


0.00029 


5.69449 


0.00023 


K (km s-l) 


6.109 


0.014 


6.113 


0.009 


(BJDtdb-2454000) 


868.9878 


0.0042 


868.9877 


0.0016 


ecoscij 


0.0000 


0.0015 


0. 




esinu) 


-0.0005 


+0.0011 

-0.0019 


0. 




e 


0.0017 


+0.0019 
-0.0017 


0. 




u! (rad) 


4.69 


1.58 


TT 

2 




mjnini (Mj) 


64.3 


3.0 


64.3 


3.0 


Systemic velocity 70 (km s"'^) 


18.67 


0.26 


18.68 


0.24 


KPNO offset (70 - Jkpno, km s-i) 


14.89 


0.26 


14.90 


0.25 


HET offset (70 - jhet, km s-l) 


12.59 


0.26 


12.61 


0.24 


MPP slope "/mpp (km day~^) 


-0.0044 


0.0065 


-0.0037 


0.0065 


KPNO slope -ykpno (km day~^) 


0.016 


0.034 


0.013 


0.033 


HET slope jhet (km s'^ day-^) 


-0.0010 


0.00035 


-0.0011 


0.00024 


Total anv Scale Factor (MPP) 


8.64 




8.49 




Total anv Scale Factor (KPNO) 


3.43 




3.37 




Total a^v Scale Factor (HET) 


0.39 




0.30 




Combined x^/dof 


1.69 




1.55 
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Table 5. Stellar host properties of TYC 2949-00557-1. 



P cir a,ni 6^ 6 r 


Va,li_iG 


U n c G rt a,i nt y 


a (J2000)'' 


101.921152 (dog) 


06:47:41.076 {HH:MM:SS) 


(5 (J2000)'^ 


42.009332 (deg) 


42:00:33.60 (DD:MM:SS) 


B 


12.846 


0.023 


V 


12.142 


0.031 




11.750 


0.039 


Ic 


11.391 


0.043 


J2MASS 


10.820 


0.022 


H2MASS 


10.474 


0.021 


K2MASS 


10.421 


0.018 


Teff (K) 


6135. 


40. 


log (9 [cms-1]) 


4.4 


0.1 


[Fe/H] 


0.32 


0.01 


vsinl (km s"-"-) 


7. 


1. 


■^primary {-^Q ) 


1.25 


0.09 


-^primary {-^0 ) 


1.15 


0.15 


d (pc) 


413 


+66 

-57 



^BVRI magnitudes are unweighted averages from Table |6] 
Coordinates taken from the Tycho 2 Catalog llH0g et al.ll2000h 



Table 6. Measured photometry for TYC 2949-00557-1 from Hereford Arizona Observatory. 



Band'' 


Flux 


Uncertainty 




(mags) 


(mags) 


B 


12.839 


0.023 


V 


12.135 


0.031 


Rc 


11.750 


0.039 


Ic 


11.390 


0.043 


9' 


12.450 


0.020 


r' 


11.937 


0.016 


i' 


11.811 


0.028 


R Sloan 


12.853 


0.023 




12.150 


0.024 


R^Sloan 


11.749 


0.029 




11.392 


0.043 



''Measurements labeled with sioan 
are converted fro m the Sloan filter 
observations using lSmith et al. 1 1I2OO2I') 
transformation equations. 
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Fig. 1. — Lomb-Scargle periodogram for the MPP data (14 epochs). A clear and highly significant 
peak at P = 5.68 days can be seen, with a power of ~ 217, which has a false alarm probability 
based on scrambling the data of < 0.01%. Powers corresponding to false alarm probabilities of 1%, 
0.1%, and 0.01% are also shown. 
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Fig. 2. — Results from the combined MPP, KPNO and HET analysis, where the MPP errors are 
scaled by a factor of 8.49, KPNO errors are scaled by a factor of 3.37 and the HET errors are scaled 
by a factor of 0.3. The eccentricity is fixed at e = 0. MPP data are in blue, KPNO data are in 
green and HET data are in red. The systemic velocity 70 has been removed. 
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TYC 2949-00557-1. Hereford Arizona Observatory 



RMS 



2010.04.15 



2009.03.16 



LU -0.070 
a: 

LU 

Li- -0.060 - 



2009.02.27 



2009.02.21 



2009.02.19 



+4.+ I****S-<i'* 



0.00519 



0.0023 



-^--j;--. ■ . ---^ ----^ f - - 



0.0038 



0.0027 



10 



11 



Fig. 3. — Photometric observations with the Hereford Arizona Observatory telescope. The horizon- 
tal axis is elapsed time each night in hours. The star is photometrically stable with an RMS of 2-4 
mmag. The increased scatter towards the end of 2009.02.21, 2009.02.27 and most of 2010.04.15 are 
due to observing at high airmass. The vertical bars are the predicted times of ingress, mid-transit 
and egress based on the RV fit in §2.41 and an assumed transit duration of 3.3 hours. The two 
mid-transit estimates are based on the two results from the RV fitting ("C" is for e = 0, "E" is 
for the case where e is left as a free parameter). The widths correspond to the uncertainties in the 
mid-transit times. 
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T3 




2 4 6 8 10 

Period 

Fig. 4.— (Top panel) KELT-N light curve for TYC 2949-00557-1. (Bottom panel) Lomb-Scargle 
periodogram of the KELT data, showing no evidence for any significant periodicities for periods 
of P = 1 — 10 days, including the period of the RV companion (vertical dashed line) and the first 
harmonic (vertical dotted line). 
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Fig. 5.— The small filled circles show the KELT-N light curve for TYC 2949-00557-1 phased 
according to the ephemeris from the joint RV fit (see §2.4p . The larger filled squares with error 
bars are binned 0.025 in phase. 
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Fig. 6. — Probability that transits of a companion are excluded at levels of Ax^ = {9, 16, 25} 
based on the analysis of the combined Hereford and KELT photometric datasets, as a function of 
the radius of the companion. The black, long-dashed line is a case where boxcar-shaped transits 
were used as a test, and is for Ax^ = 16. Transits of companions with radius r > 0.75 Rj can be 
excluded at the 95% confidence level. 
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Fig. 7. — Top: map in Tcs — Ay space showing the degeneracy between extinction and Tcff. 
Inner contours are 1-cr uncertainties for T^q and Ay assuming they are the only significant degrees 
of freedom, outer contours assume all four parameters are significant. Bottom: NextGen model 
overplotted on the observed fluxes. The dotted curve is a blackbody SED for the best-fit T^s. 
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Fig. 8. — HET/HRS template spectrum in black with best-fit model in red for the Mgb region. 
The spectra have been continuum normalized and the relative flux density is plotted against the 
wavelength in Angstroms. 




Fig. 9. — A portion of the ARCES spectrum (R ~ 31,000) used in the SME analysis. The input 
spectrum is in white with the best-fit model overlaid in black. 
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Fig. 10. — HR diagram as a function of T^s and log g based on Yonsei-Yale stellar evolution models. 
The solid track is for the best-fit stellar parameters of 1.25 Mq and [Fe/H] = +0.32. The two dashed 
tracks are for masses of 1.25 it 0.09Mq and represent the 1-a uncertainties on the mass. Blue dots 
are the location of the star at ages of {1,2,3} Gyr, respectively. TYC 2949-00557-1 is consistent 
with a ZAMS star that has an age no older than ~ 2 Gyr. 
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Fig. 11. — Cumulative probability that the mass of the companion to TYC 2949-00557-1 is less 
than a given mass in solar masses. These probabilities account for the uncertainties and covariances 
between the parameters of the radial velocity fit, the uncertainty in the mass of the primary, assum- 
ing a uniform distribution of cos i, and adopting various priors for the distribution of companion 
mass ratios dN/dlogq. 



90 67 43 20 
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Fig. 12. — Contours of the time required to synchronize the primary's rotational period to the 
orbital period in Gyr. Solid contours correspond to the best fit, dotted contours are for an unlikely 
but plausible case which maximizes the synchronization timescale. 



